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ABSTRACT 

In the past two years, tremendous progresses to understand X-ray Flashes have been 
made. Now it is widely believed that X-ray Flashes and Gamma-ray Bursts are in- 
trinsically the same, their much different peak energy and flux may be just due to 
our different viewing angles to them. Here we analytically calculate the very early 
afterglow of X-ray Flashes, i.e., the reverse shock emission powered by the outflows 
interacting with the interstellar medium. Assuming z ~ 0.3, we have shown that typi- 
cally the R Band flux of reverse shock emission can be bright to 16 ~ 17th magnitude 
(the actual values are model dependent and sensitive to the initial Lorentz factor of 
the viewed ejecta). That emission is bright enough to be detected by the telescope 
on work today such as ROTSE-III or the upcoming UVOT carried by Swift Satellite, 
planned for launching in later 2004. 

Key words: X-rays: general — Gamma-rays: bursts — radiation mechanism: non- 
thermal 



1 INTRODUCTION 

X-ray Flashes (XRFs) have received increasing attention in 
the past several years (e.g., Heise et al. 2002; Kippen et 
al. 2002). In many respects, XRFs are similar to "classical" 
Gamma-ray Bursts (GRBs): (1) Their sky distribution is 

' nearly isotropic; (2) The red-shift z = 0.251 of XRF 020903 
(Soderberg et al. 2003) and afterglows of XRF 020903 and 

■ XRF 030723 have been detected (Soderberg et al. 2002; 
Prigozhin et al. 2003), which suggests that XRFs are also 
cosmic events; (3) They have Tgo durations ranging from 
20 to 200 seconds; (4) Their spectrum can also be fitted by 
the Band spectrum (Kippen et al. 2002); (5) They have the 
temporal structure which is very similar to the X-ray coun- 
terparts of GRBs (Heise et al. 2002). However, their spectral 
peak energies -E pca k,obs ~ lOkeV are much lower than that of 
GRBs (~ 300keV). These similarities lead to the suggestion 
that XRFs are the extension of GRBs and X-ray rich GRBs 
to even softer regime (Kippen et al. 2002; Lamb, Donaghy 
& Graziani 2004). 

Now there are several models proposed to account for 
XRFs, i.e., the off-beam uniform jet model (e.g., Ioka & 
Nakamura 2001; Yamazaki, Ioka & Nakamura 2002); the 
wide opening angle uniform jet model (e.g., Lamb et al. 
2004); the Gaussian jet model (e.g., Zhang & Meszaros 
2002a; Lloyd-Ronning et al. 2004; Zhang et al. 2004a); the 
power-law jet model (e.g., Meszaros, Rees & Wijers 1998; 
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Jin & Wei 2004); the two component jet model (e.g., Zhang, 
Woosley & Heger 2004; Huang et al. 2004) and the Cannon- 
ball model (Dado, Dar & De Rugula 2003) . 

In this Letter, instead of taking a further insight into 
these models, we turn to calculate the possible accom- 
panying very early afterglows. We are mostly incited by 
the upcoming Swift Satellite^, which carries three main 
telescopes — The Burst Alert Telescope (BAT), the X-ray 
Telescope (XRT), the Ultraviolet and Optical Telescope 
(UVOT). The energy range of BAT is 15-150 keV. Con- 
sidering its high sensitivity, XRFs with not too much lower 
peak energies can be detected as well as GRBs. BAT will ob- 
serve and locate hundreds of bursts per year to better than 
4 arc minutes accuracy. Using this prompt burst location 
information, Swift can slew quickly to point on-board XRT 
and UVOT at the burst for continued afterglow studies. The 
spacecraft's 20—70 seconds time-to-target means that about 
~100 GRBs+XRFs per year (about 1/3 of the total) will be 
observed by the narrow field instruments during the gamma 
ray emission. The UVOT is sensitive to magnitude 24 in a 
1000 seconds exposure (For a linear increasing of the sen- 
sitivity with the exposure time, that means a sensitivity of 
magnitude 19 in a 10 seconds exposure). Then, though the 
very early afterglow of XRFs may be dimmer than that of 
GRBs, they can be detected directly. I 



t http:/ /swift. gsfc.nasa.gov/science/instruments 

t By now, there is only one upper limit (mpt > 19) at a time 
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2 THE VERY EARLY AFTERGLOW OF GRBS 

In the standard fireball model for Gamma-ray bursts, the 
very early afterglow of GRBs powered by the ejecta inter- 
acting with the interstellar medium (ISM) or stellar wind 
has been discussed in great detail (e.g., Sari & Piran 1999; 
Meszaros & Rees 1999; Kobayashi 2000; Li et al. 2003; Nakar 
& Piran 2004; Wu et al. 2003; Kobayashi & Zhang 2003a). 
Recently, the reverse shock (RS) emission powered by mag- 
netized outflows — medium interaction has been discussed by 
Fan, Wei & Wang (2004) and Zhang & Kobayashi (2004) 
independently. So far there are three very early afterglows 
having been detected (see Sari & Piran 1999; Kobayashi & 
Zhang 2003b; Wei 2003 and the references therein). With the 
launching of Swift, that number may be increased greatly. 

Modelling the very early afterglow can be used to con- 
strain some poorly known physical parameters, such as the 
initial Lorentz factor of the outflow, ee (the fraction of the 
internal energy converted into magnetic energy) and so on 
(e.g., Wang et al. 2000). Interestingly, it is found that the 
RS emission regions of GRB 990123 and GRB 021211 arc 
likely to be magnetized (Fan et al. 2002; Zhang, Kobayashi 
& Meszaros 2003; Kumar & Panaitescu 2003). 



3 THE VERY EARLY AFTERGLOW OF XRFS 

One of obstacles we encountered in the current work is the 
poorly known angular distribution of the initial Lorentz fac- 
tor of these jets except the off-beam one. Kumar & Gra- 
not (2003) have performed a numerical investigation on the 
hydrodynamical evolution of a Gaussian jet by assuming 
rj(9) is also Gaussian distribution. However, if the viewed 
Lorentz factor rj(9 v ) is significantly lower than 100 and 
XRFs are powered by internal shocks, the observed spec- 
trum should be thermal, which is inconsistent with the cur- 
rent observation. In their Monte Carlo simulation, Zhang 
et al. (2004) have taken r){9) as a free function but found 
that a flat (or at most slightly variable) angular distribution 
of Lorentz factor is indeed required to interpret the cur- 
rent observations, particularly, the empirical relationship^ 
Speak oc -E^ 2 • Therefore, and partly for convenience, we 
assume r/(9 v ) = const ~ 150. For the off-beam jet model, 
XRFs are intrinsically the GRBs. As usual, we take r\ — 300. 

In this work, we assume r/(9 v ) = const ~ 150 for all 
on-beam jets (i.e., the Gaussian jet, the power-law jet, the 
wide opening jet and the two component jet). Different 
from the off-beam one, for these four type jets, there are 
relativistically moving materials beaming towards the ob- 
server. As long as the Lorentz factor is large enough, within 
the I/7 cone the jet structure effect is not significant, all 
these four models would give a rather similar early afterglow 
lightcurves. Their differences would only appear in much 
later epochs when the jet-structure effect becomes promi- 
nent. So, for an illustrative purpose, in this Letter, only the 

~ 50s after the main burst of XRF 030723 has been reported 
(Smith ct al. 2003) 

§ Its validity for GRBs and XRFs has been verified by many 
authors (e.g., Lloyd-Ronning, Petrosian & Mallozzi 2000; Amati 
ct al. 2002; Zhang & Meszaros 2002b; Wei & Gao 2003; Lamb ct 
al. 2004; Liang, Dai & Wu 2004). 



possible RS emission powered by the Gaussian jet and the 
off-axis jet interacting with ISM have been calculated. 

3.1 The RS emission of the off-beam jet 

For a uniform jet, if the line of sight is slightly beyond the 
cone, i.e., A9 = 9 V — 9j ct > (#j c t is the opening angle of 
jet, 9 V is our viewing angle to the jet), the observed peak 
energy decreases as 

i^ of f = au on , (1) 

where v on / oH are the observed frequency on/off-axis respec- 
tively, a « [1 + 7 2 (A#) 2 ]~ 1 (e.g., Granot et al. 2002). On 
the contrary, the observed duration is increased by a fac- 
tor of a -1 . Then a "classical" GRB will be detected as a 
long lasting XRF (e.g., Ioka & Nakamura 2001; Yamazaki 
et al. 2002). By assuming an accompanying supernova, it is 
claimed that the off-beam jet model can fit the later after- 
glow of XRF 030723 quite well (Fynbo et al. 2004). 

The afterglow of off-beam jet has been numerical calcu- 
lated by many authors, and an empirical formula has been 
proposed to estimate the observed flux (e.g., Granot et al. 
2002; Jin & Wei 2004) 

3 

F, off (AMo bs )^yF, on (0,t). (2) 

where dt b B = dt/a. Different from the long lasting forward 
shock (FS), the RS disappears after it crosses the ejecta. 
The crossing time is estimated by t x = max{idcc, 7d ur }, 
where two timescales are involved: (i) The deceleration time 
tdec, which can be calculated as follows: The outflow is de- 
celerated significantly at the deceleration radius (Rees & 
Meszaros 1992) 

R dcc « 5.6 x 10 16 cm E^^ 3 ^, (3) 

where _Ei so ~ 10 53 ergs is the typical isotropic energy of the 
GRB outflow, m ~ 1 is the number density of the ISM, 
f] ~ 300 is the initial Lorentz factor of the outflow at the end 
of 7— ray emission phase. Throughout this Letter, we adopt 
the convention Q x = Q/10 x for expressing the physical pa- 
rameters, using cgs units. At 7?dec, the Lorentz factor of the 
outflow drops to 7d cc — ry/2. The corresponding timescale is 

£dcc ~ i?dec/27d ec C = 40 S i?doc,16.75 7dec,2. 18 ■ ( 4 ) 

(ii) 2d ur , the local duration of the GRB correspond- 
ing to the XRF, can be estimated as follows: In the off- 
beam jet model of XRFs, the observed duration of XRF is 
~ (1 + 2)a ( 7 1 T dur , where ao = [l + r/ 2 (A9) 2 ]~ 1 . Here we take 
ao — 0.03 — For much smaller ao, the observed luminosity 
(L oc o,q) is too dim to be detected; for much larger ao, the 
observed peak energy (oc ao) is in the hard X-ray energy, 
i.e., what we observed is X-ray rich burst rather than XRF. 
One potential problem of the off-beam model is that, in prin- 
ciple, the typical duration of XRFs should be tens times of 
that of typical GRBs — which has not been supported by the 
present observations (The observed duration of the XRFs 
ranging from 20 s to 200 s). Here we take T dur ~ 10s, which 
matches that of the "classical" long GRBs (~ 20s/(l + z'), 
z' ~ 1 is the typical red-shift of GRBs). 

For tdoc > T dur , the shell is thin, otherwise the shell is 
thick (Sari & Piran 1999; Kobayashi 2000). For the param- 
eters taken here, we have tdec > Tdur, so the shell is thin (If 
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the shell is thick, the following discussion is invalid and we 
refer the reader to see §3.2 for detailed treatment). Conse- 
quently, t x = idoc, 7x = 7dcc and R x = -Rdcc (7x is the bulk 
Lorentz factor of the outflow at t x , R x is the corresponding 
radius) . 

The Lorentz factor of the shocked outflow relative to 
the initial one is 



734, x « (vhx + 7xA?)/2 = 1.25, 



(5) 



which suggests that the RS is only mildly relativistic. 

At R x , all the electrons contained in the outflow have 
been heated by the RS and distribute as dn/d~/ c oc j~ p 
for 7o > 7e,m (Throughout this Letter we take p = 2.2), 
where the "minimal" thermal Lorentz factor, 7e, m , can be 
estimated by 



m p p — 2 

7c, m = —Z 7 £ c(,734,x 



m c p 



1) = 23F, 



(6) 



/ 734, X — 
-0-5 1 0.25 



■), e e is the fraction of thermal 



where F = 

energy obtained by the electrons, m p (m e ) are the rest mass 
of proton (electron) respectively. The typical synchrotron 
radiation frequency can be estimated by 



"ll,mlxeB 



2(1 + z)nm c c 



4.0 x 10 1 
1 + z 



Hz F e B,-l n l,07x,2.18, 



(7) 



where B « 0.12G nl^gf (e.g., Fan et al. 2002) 
is the magnetic strength generated in the RS. Here we 
take £b ~ 0.1 rather than 0.01 since the ejecta is likely 
to be magnetized (e.g., Fan et al. 2002; Zhang et al. 
2003). The cooling Lorentz factor of the shocked electrons 
is (Sari, Piran & Narayan 1998, hereafter SPN) -y c>x ~ 
67rm c c/(aT7xB 2 ix) ~ 360 (Assuming the involved Q x = 1), 
where <tt is the Thompson cross section. Correspondingly, 
the cooling frequency is v c>x = (^ £i2L ) 2 ^ m « 10 15 /(l + «) Hz. 
For z — 0.3 (as it is general suggested for XRFs), i/ c , x is 
larger than the observer frequency ^R, b s = 4.6 x 10 14 Hz. 
Following Wu et al. (2003; see their appendix Al for detail), 
the synchrotron self-absorption frequency of the compressed 
outflow is ^ a , x ~ 3.6 x 10 1 /(1 + z) Hz (Assuming the in- 
volved Q x = 1). So the synchrotron self-absorption effect 
can not significantly change the R band spectrum, which is 
of our interest. 

At R x , the on-beam peak flux can be estimated 
by i^, m ax(on) « (1 + z)N e m c c 2 o Tlx B/l2-KeD't « 
92.3(^) Jy £i S o,53^7x, 2.18^0 ^-i^l.It.t (SPN), 
where N e = Ei so /qm p c 2 is the total number of electrons 
contained in the outflow, Dl is the luminosity distance (we 
assume H = 70km s" 1 Mpc" 1 , ft M = 0.27, fi A = 0.73). 
With equation (2), the off-axis observed energy flux can be 
estimated by (SPN) 



"R,obs(°ff) 



4* 



/ ^RjObs \ - 
i/,max(on) \ ) 



(p-l)/2 



= 0.8mJy (9a> 



^X fm, X 
(p+5) 



1 + Z.3^p 



DT 



(8) 



where a x = [1 + (7 X A9) 2 ]' 1 . If we take a = 0.03, i.e., 
r)A9 — 5.7 (see the reasons mentioned in the paragraph be- 



low equation (4)), we have a x = 1/9, i.e., 7 X AO = 2.8. The 
observed crossing time t x ,obs ~ tx /«x ~ 360(1 + z)s. At R 
band, the magnitude rag « 17, which is bright enough to be 
detected by the telescopes on work today, or the upcoming 
UVOT. 

Here, we briefly discuss the observed very early after- 
glow light curve. In the case of thin shell, the on-axis light 
curve is well approximated by (for v m < v < v c ,x, i.e., the 
slow cooling case) 



2p 



for t < t > 



^(on) OC i t _(iip +3 )/14 i for t > t 



tdec- 



(9) 



The light curve for t > t x is presented in Meszaros & Rees 
(1999; Assuming the comoving magnetic field contained in 
the shocked outflow is freezing and taking g — 3, since for 
the parameters adopted in this letter, the FS emission is 
in fast cooling). The light curve for t < t x adopted here 
is slightly different from that of Kobayashi (2000) and Fan 
et al. (2002). Below we derive it in some detail. For the 
Newtonian RS, the bulk Lorentz factor of the ejecta 7 ~ 77, 
734 — 1 oc (ru /ni) _1 oc R (Sari & Piran 1995), where is 
the comoving number density of the outflow, R ~ 2r; 2 ct is 
the radius of the outflow. Therefore 734 — 1 oc t 2 , substitut- 
ing it into equation (7) we have v m oc t 4 . On the other hand, 
/334 oc t, which results in oc t 2 , where N' e is the number of 
shocked electrons. Substituting these relations into the ex- 
pression of -F„ max ( on ) , we have Fi,, max (on) oc t 2 . Therefore the 
observed R band light curve F VR b (on) oc t 2 (t~ 4 )~ (j '~ 1 " )/ ' 2 oc 

Assuming 7 = 77(1 — t/2t x ) for t < t x , we have 
*obs « [£-§(£-l)£]*5 For t > t x , 7 = r,(t/t x )- 3/7 /2, we 



have tabs 



B(ao-l) 



„,.. ix +[l + I(^f)(^) 6/7 ]t. With these 
relations (including eqs. (2) and (9)), we obtain one sam- 
ple light curve of the RS emission powered by the off-beam 
uniform jet-ISM interaction, which has been presented in 
figure 1 (the solid line). Naturally, for t bs < iobs.x, the in- 
creasing of the light curve is more rapidly than t 2p s . For 
tabs > tobs,x , as long as ^yA9 is not much smaller than 1, the 
factor a( p + 5 )/ 2 increases rapidly. Therefore, at early time, 
the optical emission increases, rather than decreases with 
time. At much later time, 7 AS — > 1, the factor a( p + 5 )/ 2 in- 
creases only slightly. So the light curve of RS emission drops 
as oc t~ 2 s (Figure 1, the thin solid line). 

Taking es.f = 0.01, e e ,f = e c (the subscript / repre- 
sents the forward shock), we have: At t x ,obs, F„ B oba(off) ,f ~ 
0.04mJy. For t > t x , obs , F UR obs (on),f oc ^ 1/3 (e.g., SPN). 
With eq. (2) and t ohs « ^^t x + [1 + |(i=ai)(*2L)«A]t, 
the sample early light curve of the FS emission has been 
presented in Figure 1 (The dotted line) . 



3.2 The RS emission of Gaussian jet 

For a Gaussian jet, the observed isotropic energy £i so (e v ) = 
^iso(e v =o) exp(— 6 2 /26q). Typically, the observed peak en- 
ergy of XRF is about 0.03 times that of GRBs and -Ei so (e v ) ~ 
lO -3 i?i so (0 v= o). The corresponding viewing angle is 8 V ~ 
3.7#o- Taking rj(d v ) — 150, equation (3) gives R dcc ~ 
8.8 x 10 15 cm. Correspondingly, t dcc ~ 25 s, which is much 
shorter than the typical duration of the XRFs T90 ~ 100 s, 
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Figure 1. The sample very early R-band (^R. b s = 4.6 X 10 14 Hz) 
light curves powered by the outflows interacting with the ISM. 
The solid/dotted lines represent the RS/FS emission as a func- 
tion of observer time of the off-beam uniform jet; the thick 
dashed/dash-dotted lines represent the RS/FS emission as a func- 
tion of time of the Gaussian jet. The parameters for plotting the 
solid/dotted lines arc: T dur = 10s, A6> = 0.019rad, r) = 300, 
E iso = 10 53 crgs, z = 0.3, p = 2.2, e c = e e ,f = 0.3, e B = 0.1, 
e B,f = 0.01. For plotting the thick dashed/dash-dottcd lines, 
the parameters are the same except T90 = 100s, r](8 v ) = 150, 



,) = 10 bu ergs. 



i.e., the shell is thick. So, locally t x w Tgo/(l + z), with 
which R x and y x can be calculated self-consistently. 

At R x , the energy conservation of the system, i.e., the 
shocked ISM and the shocked viewing outflow, gives 



7x734, x M cj + 7 X M SW « jj(0 v )M ej , 



(10) 



where M c j (M sw ) is the mass of the viewed ejecta (the swept 
ISM). In the thick shell case, the RS is mild-relativistic 
and 734, x ~ 7?(#v)/27x- Now equation (10) reduces to 
7xM sw « r?(6> v )M cj /2. Considering that M sw = ^i? 3 nim p 
and Tgo/(l + z) ~ i? x /27 x c, we have 



R> 



7x 



1.4 x 10 16 cm El/ o 4 iw n lfi 



1/V/ 4 



'90,2 



1/4 



55-Ey^,, , rn n 

iso(y v ) ,50 



-1/8^-3/8/ 1+£n 3/8 
1,0 J 90,2 I y J 



(11) 



(12) 



Now 734, x ~ 1.55, which is mild-relativistic. So the assump- 
tion made before is reasonable. Similar to section 3.1, the 
typical frequency of RS emission can be estimated by 



^m, x 



3.8 x 10 1 



1 + z 



Hz i ? ini 2 ne B ,-l7x,1.74- 



(13) 



where F\ 
and z — 



i - 

e e ,-o.5{- 



0.55 

0.3 we have v c 



Similarly, taking Q x = 1 
- 2.9 x 10 16 Hz. Follow- 



ing Wu et al. (2003), the synchrotron self-absorption fre- 
quency is i/ a ,x ~ 4.5 x 10 10 Hz. Therefore both of them 
can not affect the R band spectrum significantly. The 
peak flux of RS emission can be estimated by i^, ma x ~ 



28mJy {^)E iso{ g v) ^ r) 2 i sl l^. 7A n\^ 4(-i£>l,2 
the observed peak energy flux can be estimated by (SPN) 



Then 



^R.obs ^-(p-l)/2 



= 1.2mJy S iso(£)v)i50 772.i87x,i.74 i? i P e B ,-i n i,o 

(^J^^.r, (14) 

the magnitude m,R w 16, which is also bright enough to be 
detected by the upcoming UVOT or the telescopes on work 
today, as long as the response to the XRF is fast enough. 

In the case of thick shell, the very early light curve of a 
standard fireball takes the form (For the parameters adopted 
in this Letter, the reverse/FS emission are all in slow cooling, 
for t > t x ). 



\ t - 3 < 



3(5p+l)/16 



for t < tx 
for t > tx 



■T 90 /(l + z), 
■T 90 /(l + z). 



(15) 



Thanks to the beaming effect, these scaling laws may 
be applied to the current work as well. Here we simply 
take equation (15) to plot the sample very early light curve 
powered by Gaussian jet — ISM interaction (See figure 1, 
the thick dash line). The corresponding light curve of the 
FS emission (see figure 1, the thick dash-dotted line) is 
plotted as follows: Following SPN, at t x ,obs = (1 + z)t x , 
F -R,obs,f ~ 0.09mJy. For t < t , F„ R oba ' f tx i 1 / 2 (where 
to is determined by f m ,f(to) = (1 + z)u R . obs ); For t > to, 

As shown in Figure 1, in many respects, i.e., the peak 
time and the temporal behavior, the current light curve is 
much different from that powered by the off-beam jet — ISM 
interaction, which may help us to distinguish them. For ex- 
ample, the FS peak emission of the off-beam jet is much 
brighter than that of Gaussian jet, which is mainly due to: 
In the off-beam jet model, we take A8 — 0.019rad. At several 
thousand seconds after the main burst, the ejecta has been 
decelerated significantly (7 < 50), as a result, the beaming 
effect is not important. So the viewed isotropic energy of 
the "off-beam" ejecta is nearly Ei BO , which is much larger 



than that of the Gaussian jet model (~ 



so(9 v )J 



Tn which 



9 V ~ 0.2rad (Zhang et al. 2004), the emission powered by 
the central energetic ejecta can not be observed at early time 
as long as 7 ^> 5. 



4 DISCUSSION AND CONCLUSION 

In the past several years, XRFs have received more and more 
attentions, but their nature remains unknown. Considering 
the similarity between XRFs and GRBs on duration, tempo- 
ral structure, spectrum and so on, these bursts may be the 
same phenomenon. The different peak energy as well as the 
peak energy flux may be just due to our different viewing 
angles to them. This viewpoint has been supported by the 
recent detection of two afterglows and one red-shift of XRFs. 
However, there are more than 6 models have been proposed 
to explain the XRFs. Some of them (for example, the Gaus- 
sian jet model, the off-beam uniform jet model) work well 
on explaining the current observation. Perhaps only the de- 
tailed multi- wavelength afterglows modelling (including the 
very early afterglow discussed in this Letter) can provide us 
a reliable identification on them. 

With two current leading models, in this Letter, the 
very early RS emission of XRFs have been analytically in- 
vestigated. As XRFs are much dimmer than GRBs, the pre- 
dicted very early R band afterglow of XRFs is dimmer than 
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those of GRBs. But some of them, if not all, are still bright 
enough to be detected by ROTSE-III on work today or by 
the upcoming Swift mission. The actual results are model 
dependent, which may in turn provide us a chance to see 
which one is better, if the very early afterglow of XRFs has 
been really detected. One thing should be emphasized here 
is that the predicted R band flux is sensitive to the initial 
Lorentz factor, i.e., -F„ R oba oc 7?(0v) _1 7x +1 (734, x - l) p ~\ 
So, if r)(9 v ) ~ tens rather than 150 taken in this Letter, the 
resulting very early optical emission will be much dimmer. 

As realized by more and more "GRB people" , modelling 
the very early afterglow of GRBs can impose some strin- 
gent constraint on the fundamental physical parameters of 
the outflow such as the initial Lorentz factor of the ejecta 
(e.g., Sari & Piran 1999; Wang et al. 2000), or help to see 
whether the outflows are magnetized or not (e.g., Fan et al. 
2002; Zhang et al. 2003; Kumar & Panaitescu 2003; Fan et 
al. 2004b; Zhang & Kobayashi 2004). In addition, the RS 
emission powered by the ejecta — stellar wind interaction is 
much different from that powered by the ejecta — ISM in- 
teraction (e.g., Wu et al. 2003; Kobayashi & Zhang 2003a; 
Fan et al. 2004b). Therefore the very early afterglow obser- 
vation can provide us an independent chance to determine 
the environment where the "classical" GRB was born in. In 
principle, interstellar wind environment may be common, 
and the predicted very early R band emission is very strong 
(ma ~ 9 or even brighter). So far there are only three very 
early afterglow having been reported. All of them can be 
well fitted by the ejecta — ISM interaction model (e.g., Sari 
& Piran 1999; Kobayashi & Zhang 2003b; Wei 2003). There- 
fore, in this Letter only the ejecta — ISM interaction case has 
been investigated. Fortunately, it is straightforward to ex- 
tend our treatment to the wind case. Anyway, the impor- 
tance of modelling the very early afterglow of GRBs can be 
applied to that of XRFs as well. 

In this Letter, the possible pair loading during the 7/X- 
ray burst phase has not been taken into account. At least 
for the off-beam uniform jet model, during the initial 7— ray 
emission phase, large amount of electrons/positrons may be 
created. The annihilation time scale is long and most of gen- 
erated pairs can not be annihilated locally. These pairs will 
be heated by the RS, too (Li et al. 2003, and the references 
therein). Without doubt, the RS emission has been signif- 
icantly softened. But, as argued by Fan, Dai & Lu (2004), 
the R band emission may be not. This can be easily un- 
derstood, if there are k times electrons/positrons of that 
associated with the baryons, current v m would be v m /k . 
On the other hand, current F„ ima x would be fcF„, m ax. Then, 
in the slow cooling case, current F^obs should be k 2 ~ p F v , ohs . 
For p ~ 2.2, such dependence is far from sensitive. In other 
models for XRFs, during the X-ray emission phase, at least 
in the viewed area, the possible pair loading process is unim- 
portant and can be ignored safely. 
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